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Quasi-single field inflation predicts a peculiar momentum dependence in the squeezed limit of the 
primordial bispectrum which smoothly interpolates between the local and equilateral models. This 
dependence is directly related to the mass of the isocurvatons in the theory which is determined 
by the supersymmetry. Therefore, in the event of detection of a non-zero primordial bispectrum, 
additional constraints on the parameter controlling the momentum-dependence in the squeezed limit 
becomes an important question. We explore the effects of these non-Gaussian initial conditions on 
large-scale structure and the cosmic microwave background, with particular attention to the galaxy 
power spectrum at large scales and scale-dependence corrections to galaxy bias. We determine the 
simultaneous constraints on the two parameters describing the QSF bispectrum that we can expect 
from upcoming large-scale structure and cosmic microwave background observations. We find that 
for relatively large values of the non-Gaussian amplitude parameters, but still well within current 
uncertainties, galaxy power spectrum measurements will be able to distinguish the QSF scenario 
from the predictions of the local model. A CMB likelihood analysis, as well as Fisher matrix analysis, 
shows that there is also a range of parameter values for which Planck data may be able distinguish 
between QSF models and the related local and equilateral shapes. Given the different observational 
weightings of the CMB and LSS results, degeneracies can be significantly reduced in a joint analysis. 
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I. INTRODUCTION 

Quasi-single field (QSF) inflation models [1-3] are a natural consequence of inflation model-building in string theory 
and supergravity. In addition to the inflaton field, these models have extra fields with masses of order the Hubble 
parameter. Such masses are stabilized by the supersymmetry. A distinctive observational signature of these massive 
fields is a one-parameter family of large non-Gaussianities whose squeezed limits interpolate between the local and 
the equilateral shape. Therefore, by measuring the precise momentum-dependence of the squeezed configurations in 
the non-Gaussianities, in principle, we are directly measuring the parameters of the theory naturally determined by 
the fundamental principle of supersymmetry. 

The possibility of detecting a non-Gaussian component in the initial conditions of Early Universe has been the 
subject of considerable attention in recent years both from an observational perspective and theoretically through 
inflation model-building [4-6] . Current constraints from measurements of the bispectrum of temperature fluctuations 
in the cosmic microwave background (CMB) from the WMAP satellite [7] are still consistent with Gaussianity. 
However, the Planck mission [8] will soon significantly improve the errors on non-Gaussian parameters, leading to 
strong new constraints or what could be a major breakthrough in cosmology. 

The effects of non-Gaussian initial conditions on the large-scale matter and galaxy distributions have been the 
subject of several studies for more than a decade (see [11] and [5] for recent reviews). The most direct of such 
effect consists in the additional contribution to the matter bispectrum due to the linearly evolved initial component. 
Measurements of the galaxy bispectrum in upcoming, large- volume galaxy surveys are expected to improve even over 
ideal CMB limits, for any non-Gaussian model [12, 13]. 
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In addition, the relatively recent discovery of a significant scale-dependent corrections to galaxy bias due to non- 
Gaussian initial conditions [9] has led to constraints on the local non-Gaussian parameter /^l' from current Large- 
Scale Structure (LSS) data-sets, already comparable to those from the CMB [10]. Such bias corrections are present 
for models where the curvature bispectrum takes large values in the squeezed limit and precisely for this reason, the 
case of Quasi-Single Field inflation is of particular interest. The results of rcf. [9] motivated a significant number of 
further works aimed at a rigorous theoretical description of the effect [10, 14-29]. At the same time several groups 
ran new sets of simulations with non-Gaussian initial conditions of the local [18, 20, 22, 30] and other types of models, 
such as equilateral, orthogonal and folded [29, 31-34], as well as generalized local models (e.g. cubic g^L or scale- 
dependent /nl models) [24, 35-37]. The picture emerging from such extensive investigations is that the relatively 
simple expression sufficient to describe the effect of local non-Gaussianity, requires additional corrections to accurately 
describe models presenting a squeezed limit of the curvature bispectrum different from the local one [27, 29]. More 
generally for such models, the overall correction to bias induced by a generic non-Gaussian model will depend on the 
halo mass in a non-trivial way and will include a scale- independent component [18]. The extension of these predictions 
to nonlinear bias and to the description of the galaxy bispectrum have been studied in [21, 38-40] and recently tested 
in simulations by [22, 41] in the context of the local model. 

The detectability of the bias correction, both in galaxy and cluster surveys, has also been the subject of several 
works in recent years [10, 42-51]. Typical forecasted errors for upcoming galaxy redshift surveys on /^', marginalized 
over cosmological parameters, are of the order of A/^l' ~ few [46, 50], while cluster surveys including information 
on the cluster spatial correlations should provide similar results [47, 48, 51]. Beyond local non-Gaussianity, ref. [50] 
provided forecasts for equilateral and orthogonal models as well, combining information from galaxy and weak lcnsing 
observables. When only the 3D galaxy power spectrum is considered, therefore relying mostly on corrections to galaxy 
bias, the expected errors are of the order of A/^ ~ 12 and A/^ — 37 for orthogonal and equilateral non-Gaussianity, 
respectively, including CMB priors from the Planck mission. The case of an initial bispectrum described in terms 
of two parameters, and the possibility of determining them both, a case complimentary to the one considered here, 
has been considered by [52] in the context of running non-Gaussianity. In this phenomenological model, theoretically 
motivated in [53, 54], the non-Gaussian parameter / NL presents a scale-dependence parametrized in terms of a running 
parameter jing = dlnf^/dlnk. 

Along with these predictions, as already mentioned, other papers derived constraints on non-Gaussian parameters 
from current observations. In particular, ref. [10] find —29 < < 70 at 95% C.L. combining different data- 

sets in the Sloan Digital Sky Survey (SDSS), with a dominating contribution from the photometric quasar sample 
[55]. Ref. [56] finds instead the 2-cr limits 3 < /^l' < 103 and 10 < < 106 from high- redshift radio sources 

from the NRAO VLA Sky Survey (NVSS) [57] and the SDSS quasar sample [58]. More recently ref. [59] considers the 
analysis of high-redshift probes for the equilateral and folded models, in addition to the local one, finding, respectively 
-480 < /£l < 580, -7 < < 373 and 8 < f$£- < 88 at 95% C.L. It should be remarked that the model assumed 
to describe the bias correction in [59] is quite approximate as it neglects the dependence on the halo mass and further 
corrections considered for instance in [27-29], so the limits derived in the equilateral and folded cases are to be 
considered, as pointed-out by the authors, as limits on "effective" non-Gaussian parameters and /^l- 

In the first part of this work, we study the effects of QSF models of inflation on Large-Scale Structure with special 
attention given to the scale-dependent correction to the linear halo bias. The interesting aspect of this correction is 
its direct dependence on squeezed configurations of the curvature bispectrum. Since, as we will see, the momentum- 
dependence of these configuration is directly related to one important parameter of the theory, this model constitutes 
a veritable case study for non-Gaussian effects on halo bias. We will consider the relative importance of both 
scale-dependent and scale-independent corrections as a function of the halo mass. We will perform a Fisher matrix 
analysis to assess the detectability of the effect from measurements of the galaxy power spectrum in large-volume 
redshift surveys and provide a first estimate of the expected uncertainty on the two parameters controlling the initial 
bispectrum. In addition, we will consider as well the constraints on these parameters expected from measurements of 
the CMB bispectrum, with particular reference to the upcoming Planck satellite. 

This paper is organized as follows. In Section II we introduce the Quasi-Single Field model of inflation and present 
the template assumed to describe the predicted curvature bispectrum. In Section III we study the effect of the 
QSF model on the matter bispectrum and skewness, deriving some basic results useful for the following analysis. 
In Section IV we consider non-Gaussian corrections to the linear halo bias and show the results for a Fisher matrix 
analysis of the galaxy power spectrum. In Section V we discuss instead the expected constrains from CMB observations 
and discuss their combination with LSS forecasts. Finally, we present our conclusions in Section VI. 
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II. QUASI-SINGLE FIELD INFLATION 
A. Theory 

Inflation model building in supergravity and string theory naturally leads to models of quasi-single field inflation. In 
this class of multifield models, there is one field direction, the inflaton direction, that satisfies the slow-roll conditions 
through either symmetry or fine-tuning, and many other directions, the isocurvaton directions, that have masses of 
order the Hubble parameter, H. Supersymmetry plays an essential role in determining the masses of these isocurva- 
tons. Without supersymmetry, at the tree level, the coupling of the scalar fields to the space-time curvature would 
also lead to masses of order H , but these masses will run away due to loop corrections, analogous to the situation 
of the Higgs mass in particle physics. In cosmology, supersymmetry provides the only dynamical mechanism for the 
masses to maintain this value. The radiative corrections to the mass from the scalar and fermion loops automatically 
cancel down to the supersymmetry breaking scale H. Therefore, despite the Hubble parameter being determined 
by a sector independent of the isocurvatons, and no matter how large or small H is, the mass of light scalars will 
always trace the value of H through the universal gravitational coupling. For the inflaton, this mass is the origin 
of the ^-problem in supergravity inflation models, and needs to be tuned away. For the isocurvatons, they become 
signatures of supersymmetry in the primordial universe. Finding observational evidence of such scalars constitutes 
an outstanding theoretical and experimental challenge. 

The generic couplings between these isocurvatons and the inflaton have several possible consequences on the pri- 
mordial density perturbation. The simplest one is the correction to the two-point correlation function, the power 
spectrum, of the density perturbation. Since the power spectrum is a function of one momentum, the observable 
effects only appear if such a correction is non-scale-invariant. Furthermore, we expect these effects easily can be de- 
generate with other types of corrections. The most distinctive signatures of these isocurvatons come from the three or 
higher-point correlation functions. Unlike the inflaton, the self-interactions of these isocurvaton fields are free of any 
slow-roll conditions and can be very large. These become the sources of large non-Gaussianities. More importantly, 
these non-Gaussianities turn out to have very special properties. For example for the scalar three-point function ($ 3 ) 
in the simplest QSF model, the momentum dependence in the squeezed limit £3 <C k\ = fc 2 is given by [1, 2], 



($3) 



k- 2+a 0< — < 9 + (\n^ 



3/ 2 !„ «3 



fc 3 m 2 9 (1) 



fci ' H 2 ~ 4 ' 



where 



1 9 m 2 

This is a property of the shape of the non-Gaussianity and is present even if the non-Gaussianity is perfectly scale- 
invariant. For QSF, the stability of the inflaton requires m 2 > 0, so that a > — 1. As m 2 / H 2 > 9/4, the isocurvatons 
gradually become too massive to have significant effects on density perturbations in the absence of sharp features. 
So we are mainly interested in —1 < a < 1/2, i.e. < v < 3/2. Such a momentum dependence lies between that 
of the equilateral shape (a = 1) which arises in single field models with large non-Gaussianity or its direct multifield 
generalization, and that of the local shape (a = —1) which arises in multifield models with light isocurvatons m <C H 
(see [6, 60, 61] for reviews.) 

While the detailed dependence of a on the isocurvaton masses may be model-dependent in more general situations, 
the signature intermediate momentum dependence in the squeezed limit is a robust evidence for the existence of such 
isocurvatons. This can be seen qualitatively as follows [1, 2]. The fluctuations of the massive scalars decay after 
the horizon-exit. For heavy scalars they decay immediately after the horizon-exit, and for lighter scalars they decay 
more slowly. The scalar interactions, responsible for the large non-Gaussianities, are therefore generated between the 
horizon scale and the superhorizon scales. The former is responsible for the equilateral-like shapes, and the latter 
for local-like shapes. As a consistency check, if we look at the special limit of massless scalars, the superhorizon 
fluctuations do not decay, and we recover the characteristic local shape in the squeezed limit. This momentum 
dependence can be also seen more quantitatively as follows [3] , at least for a close to — 1 . Ignoring the physics within 
and near the horizon scale, the squeezed limit of the three-point function can be regarded as the modulation of the 
two-point function of two short-wavelength modes from a long-wavelength mode. After horizon exit, we know that 
the amplitude of a massive scalar decays as ~ a^ 1 ^" as a function of the scale factor a. So the amplitude of the 
long- wavelength mode has decayed by a factor of (ks/ki) 1+a by the time the short-wavelength modes start to exit the 
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horizon. The amplitude of the modulation is proportional to the amplitude of the long-wavelength mode. Taking the 
massless limit as the reference point, at which we know the squeezed limit behavior ~ fc^ 3 from the simple locality 
argument [62-64], in the massive case we get the momentum dependence <~ k^ 2+a . 

To illustrate how close these intermediate shapes can get to the local shape but with qualitatively different values 
of the fundamental parameter, we look at the example of m = 0.5 H. This mass is still of order H, qualitatively 
different from the massless isocurvaton (m <C H) in multifield slow-roll inflation models. But the resulting shape has 
the momentum dependence ~ k^ 29 , very close to the local shape k^ 3 characteristic of the massless isocurvatons. 
Therefore, how well we can measure the squeezed limit behavior, i.e. the parameter v, is an important question, if 
any large non-Gaussianities are discovered. 



B. Initial curvature bispectrum 



We will assume throughout the following template for the bispectrum of the Bardeen potential $ (with $ = 3£/5) 

[2] 

B^(k 1 ,k 2 ,k 3 ) = 6ClF(k 1 ,k 2 ,k 3 ), (3) 

with 

F(fcl ' fe ' fc3) - AW7) /NL ( fclW /2 (fcl + fc2 + fc3) 3/ 2 • W 

where N v is the Neumann function of order v. For simplicity, we will further assume scale invariance for the curvature 
power spectrum 1 , with the constant C$ defined as P$(fc) = C$/fc 3 . From now on we will generically refer with /nl to 
the non-Gaussian amplitude parameter for QSF models as described by the template above. If need be, non-Gaussian 
parameters for other models will be denoted explicitly with a superscript as, for instance, /^l' or for the local 
and equilateral shapes, respectively. 

As shown by [2], the template of Eq. (4) well reproduces the main features of the bispectrum predicted by QSF 
models, and, in particular, provides the correct scale-dependence for squeezed configurations as a function of the 
parameter v. In fact, in the squeezed limit, i.e. fc 3 <C fci, fc 2 , the leading order expression for the template becomes, 
for v £ 0, 

r ik k k \ k3<< k' k2 f r 2 18 (fci + fc 2 ) 3 - 3/2 l , . 

B^k u k 2 M) - -f™ c * A „ nNi/{8/27) (fclfc2 )3/ 2+ , k m + „- (5) 

For the special case v = 0, the limit is given by 

D „ , I ,fe 3 «fci 1 fe2 , ^ 2 36\/3 7 + ln[4fc 1 fc 2 /(fc 1 +fc 2 ) 3 ] + In fc 3 1 

B,( kl ,k 2 ,k s ) - -fa C *^Y6W0 (fcifc 2 ) 3 / 2 (fci + fc 2 ) 3 / 2 jg7»> ^ 

7 being the Euler constant 2 . In particular assuming fci = k 2 = k s and fc 3 = k, we have for v ^ 0, 



fc<<fe 3 t „ 2 9V3 2T(i/) ks 



-9/2 



while for v = we have 



S $ (fc s , fc s , fc) -/nlCI ^^ /27) ( 7 - In 2 + In k) ^ , (8) 

In the upper panels of Fig. 1 we compare the squeezed limit of flattened triangular configurations of the curvature 
bispectrum, that is B$(k, k s — k, k s ) with k s constant as k — > 0, as described by the template above with the numerical 



1 Our Fourier transform convention implies the following definitions for power spectra and bispectra {^kx^k^) = <5zj(kl2) P$(fci) and 
<*ki*k 2 *k 3 > = S D (k 12 3) B i ,(k 1 ,k2,k 3 ) with k 4j = + k,-. 

2 For v < 0.5 the next to leading order term goes like ~ fc 3 / 2_I/ , therefore for u = should be taken into account. 
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FIG. 1: Upper panels: comparison of the bispectrum B<s> as described by the template, Eq. (4), with the numerical evaluation 
for flattened triangles, B(k s , k, k s — k) as a function of k with constant k s — 1 /iMpc" 1 , for v = 1, 0.5, 0.3 and (left to right). 
Lower panels: same comparison for squeezed isosceles triangles, B(k,k s ,k s ) as a function of k with k s — 1/iMpc -1 . Notice 
that a lesser number of points has been evaluated in the v = 1 case. 



evaluation performed in [2] for v = 1, 0.5, 0.3 and 0. The lower panels of Fig. 1 show the same comparison but for 
squeezed, isosceles triangles, i.e. B^(k, k s , k s ). It is evident that the template provides an accurate description of the 
numerical results for v = while it presents an increasing (but asymptotically constant) discrepancy for larger values 
of v, of fews tens of percent in the squeezed limit. We remark that this is a minor problem as long as the templates 
provides the corrected momentum dependence in the squeezed limit. As discussed, this momentum dependence is 
the most general consequence of the models and is related to the most interesting underlying fundamental physics. 
In the squeezed limit, this discrepancy simply rescales the definition of /nl- Nonetheless, more accurate comparisons 
between observations and a specific theoretical model will clearly require a direct evaluation of the bispectrum from 
the model. So far for the simplest models, direct analytical expressions for the entire shapes are not written in a 
closed form, and explicit evaluation is only available numerically. 

A "local limit" is attained for v — 3/2, leading to B^(k^ -C k\, kz) ~ k^ 2 , as for the usual local model. It should 
be noted, however, that even for v = 3/2 and for squeezed configurations, the amplitude of the QSF template docs 
not coincide with the amplitude of the local model for the same configurations. The upper panels of Fig. 2 show 
a comparison between the QSF template for v — 3/2 and the local model bispectrum for both flattened (left) and 
isosceles (right) squeezed configurations. The shape function F(ki,k2,ka) is normalised in such a way to obtain 
F(k,k,k)k 6 = 1 for any value of k, as it happens for the local model, hence the agreement for k = 1/iMpc -1 in 
the plots. For all other configurations this is generically not true. Even for the particular case of v = 3/2 we will 
denote by /nl the amplitude parameter for the QSF bispectrum, while we will indicate with f^£ the analogous 
parameter for the local model. We find that for squeezed configurations, the two model present the same amplitude if 
/nl — 2.35/jJl'j with the factor of 2.35 representing the discrepancy shown in the plots as k — > 0. This relation is useful 
in comparing the errors on the parameters expected from observations of the galaxy power spectrum since, as we will 
see in section IV, the scale-dependent non-Gaussian corrections depend almost exclusively on squeezed configurations. 
For instance, an expected error of the order of A/^l' = 5 would correspond to an error on /nl of A/nl — 12. More 
importantly, current LSS constraints —29 < f^£- < 70 at 95% C.L. [10] translate to —68 < /nl < 164 for v = 1.5, so 
that values of /nl = 150 are still allowed by galaxy power observations within the 2-er limit. Comparisons with CMB 
results are less simple, as we will see, since the CMB bispectrum is naturally sensitive to all triangular configurations. 

The lower panels of Fig. 2 compare instead the QSF template for v — 1/2 with the "orthogonal" template proposed 
by [65] in their Eq. (3.2). It should be noted that such template represents a good approximation to the exact 
orthogonal shape for configurations far from the squeezed limit, and it is therefore a viable choice only for CMB 
analysis. Its peculiar squeezed limit, in fact, leads to a 1/k scale-dependent correction to galaxy bias, not predicted 
by the true orthogonal shape, but of some phenomenological interest, as it determines an effect on bias intermediate 
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FIG. 2: Upper panels: comparison of the QSF template for v = 1.5 with the local shape for flattened squeezed, B(k s , k, k s — k) 
(left panel) and isosceles, squeezed configurations, B(k,k s ,k s ) (right panel) assuming k s = 1/iMpc -1 . Lower panels: comparison 
of the QSF template for v — 0.5 with the orthogonal shape for the same triangular configurations. 



between the local and equilateral models [29, 32]. On the other hand, QSF inflation naturally predicts this kind of 
behavior. As shown by the figure, for v = 0.5 the asymptotic behavior in the squeezed limit for the two models is 
the same, with a difference in amplitude of about a factor of two. Notice that for the orthogonal case we assume 
the same absolute value for /nl as the QSF model, but opposite sign as the orthogonal bispectrum is negative for 
such triangles. Far from the squeezed limit, as it is evident in particular from the plot showing isosceles triangles, 
the two templates are significantly different. For this reason, we will not consider further comparisons between the 
approximate orthogonal template, but we will instead confront QSF predictions with those of both the local and 
equilateral shapes. 

For all the subsequent calculations in this paper we will assume a flat ACDM cosmological model with the following 
parameters: h = 0.719, ATI 2 = 2.41 x 10~ 9 , il b = 0.0441, il c = 0.214 and n s = 1, leading to a 8 = 0.86. The matter 
transfer function is computed with the CAMB code 3 . Notice that, for simplicity we are assuming scale invariance, 
consistent with the original expression for the bispectrum template proposed by [2]. 



III. THE MATTER BISPECTRUM AND SKEWNESS 



The most direct effect of primordial non-Gaussianity on large-scale structure is given by its linear contribution to 
the matter bispectrum. At large-scales, in fact, it is possible to approximately describe the matter bispectrum by 
its tree-level expression in Eulerian Perturbation Theory (see [66] for a general review on perturbation theory and 
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Matter bispectrum, B(k, k s , k s ), k s = 1 h Mpc 
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FIG. 3: Left panel: comparison of the primordial component to the matter bispectrum due to QSF models (red, dashed curves) 
to the component due to the local (black, dot-dashed curve) and equilateral (black, dotted curve) non-Gaussian models as well 
as to the one due to gravitational instability (gray, continuous curve) for squeezed triangular configurations, i.e. B(k,k s ,k a ) 
as function of k for fixed k s = l/iMpc^ 1 . We consider the cases given by v = 0, 0.5, 1 and 1.5 corresponding to the 
increasingly long-dashed curves from bottom to top. Right panel: similar comparison for the reduced matter bispectrum with 
fixed ki = 0.01 /iMpc -1 and k 2 = 1.1/iMpc -1 as a function of the angle 9 between ki and k 2 : now the different curves 
correspond to the same models, but including the gravity contribution. 



[5] for the case of non-Gaussian initial conditions). This is given by the sum of the primordial component and the 
contribution induced by gravitational instability. We have 

B(k 1 ,k 2 ,k 3 ) =B (k 1 ,k 2 ,k 3 ) + B G (k u k 2 ,k 3 ), (9) 

where the primordial component, linearly evolved to redshift z, is given by 

B Q (k x ,k 2 M)=M(k 1 )M(k 2 )M(k 3 )B^(kxMM), (10) 

with the function 

(id 

expressing the Poisson equation as 6^ = M(k) and 5^ representing the matter density contrast in Fourier space. 
The gravity-induced component is (see e.g. [66]) 

B G (ki, k 2 , k 3 ) = F 2 (k 1 ,k 2 )P (k 1 )P (k 2 ) + 2 perm. , (12) 

where Po(k) is the linear matter power spectrum and the second-order kernel F 2 of the perturbative expansion of 5^ 
is given by 

5 l k r k 2 (k x k 2 \ 2 / k^vkA 2 
F 2(ki,k 2 ) = - + -—— - + - +- — — . (13) 



7 2 k\ k 2 \ k 2 k\) 7 \ k\ k 

On the left panel of Fig. 3 we compare the primordial component to the matter bispectrum due to QSF models 
to the component due to the local and equilateral non-Gaussian models as well as to the one due to gravitational 
instability for squeezed triangular configurations, i.e. B(k,k s ,k s ) as function of k for fixed k s = 1/iMpc -1 . In 
particular we consider for the QSF models the cases given by v = 0, 0.5, 1 and 1.5 corresponding to the increasingly 
long-dashed curves from bottom to top. We notice how as v increases the curves approach the local model, reaching 
the same dependence on the scale for v = 1.5. The right panel shows the reduced matter bispectrum, defined as 
Q(ki,k 2 , fc 3 ) = B(k\, k 2 , k 3 )/[P(ki)P(k 2 ) + 2 perm.] where the different curves correspond to the same, different NG 
models, but where they include in all cases the gravity contribution. For the choice of triangles considered in the limit 
9—^-0, with two sides fixed at close values k\ = 0.01 /iMpc -1 and k 2 = 1.1 /iMpc" 1 . this corresponds to the squeezed 
limit, where different values of v lead to distinctly different behaviors. 
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FIG. 4: Reduced skewness S3, defined in Eq. (14), as a function of the halo mass for local (black, dot-dashed curve), equilateral 
(black, dotted curve) and QSF non-Gaussianity with v — 0, 0.5, 1 and 1.5 (red, increasingly long-dashed curves). Assumes 
/nl = 100, for all models. 



A quantity directly related to the linear matter bispectrum, and useful for the following calculations is given by the 
reduced skewness of the matter density field S R smoothed on the scale R defined as 

where the third order moment (S R ) is given by the integral 

(S%) = ( d 3 q 1 d 3 q 2 W R (q 1 )W R (q 2 )W R (q 12 )B Q (q ll q 2 ,q 12 ), (15) 



with W R (x] being the Fourier transform of a top-hat function. We notice that the reduced skewness S3 does not 
depend on redshift as well as on the normalization of primordial fluctuations and it is generically mildly depends on 
the other cosmological parameters. 4 In addition, it is almost constant even with respect to same variable R. This 
evident from Fig. 4 where we plot s$(R) for local, equilateral and QSF models with v = 0, 0.5, 1 and 1.5 assuming in 
all cases /nl = 100. 

The effect of primordial non-Gaussianity on matter correlators is not limited to the large-scale bispectrum, but 
also involves corrections to the small-scales nonlinear evolution of both power spectrum [16, 31, 67] and bispectrum 
[68, 69]. In the case of the bispectrum, such small scale corrections can be significant, but can be directly accessible 
only via weak lensing measurements. A proper assessment of the possibility offered by future weak lensing surveys 
to constrain primordial non-Gaussianity with the measurement of the shear higher-order correlation functions is not 
yet available. Instead, studies of non-Gaussian effects on galaxy correlators have witnessed recently a great deal of 
activity, mainly due to the scale-dependent corrections to galaxy bias which we will consider in the next section. 



IV. LINEAR HALO BIAS 



Only relatively recently, N-body simulations with local non-Gaussian initial conditions have shown that the bias of 
dark matter halos receives a significant scale-dependent correction at large scales [9, 18, 20, 29, 30, 32, 35]. Several 
papers assumed different approaches in the theoretical description of this effect, mostly based on the peak-background 
split framework or on the theory of peak correlations (see for instance [70] and references therein). 

In this work, we will consider the approach of [29], based on the peak-background split argument [71, 72], tested in 
N-body simulations assuming different non-Gaussian models for the initial conditions. We assume that the Eulerian 



4 In the subsequent calculation we make use of a fit for the reduced skewness of QSF models, function of the parameter v and of 
x = S c /a R given by s 3 (x,u) = exp [ci(v) lnrr + c 2 (u)(\nx) 2 + c 3 (u)(lnx) 3 + c 4 (u)(\nx) 4 } , with c x (v) = -8.61 - 0.164 v + 0.240 u 2 , 
c 2 (u) = -0.0173 + 0.0419 u - 0.0760 v 2 , c 3 (u) = -0.0351 - 0.0176v + 0.0298 u 2 and c 4 (v) = 0.00511 + 0.00330 v - 0.00513i/ 2 , accurate 
at the 1% level for 0.5 < x < 14. 
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relation, in Fourier space, between the halo density contrast 8h(m) for a given halo mass to and the matter density 
contrast 8 is given, at the linear level, by 

S h (m,k) ~ b h (m, k)5k, (16) 

with the linear halo bias function 

bh(m, k) = b si (m) + Ab sd (m, k) , (17) 

where we distinguish a scale-independent contribution b s i(m) from the scale-dependent correction Ab s d(m,k). 

In this section we will present theoretical results regarding the linear halo bias, keeping in mind that the effects 
described naturally translate into effects on the linear bias of the galaxy distribution. Galaxy bias can in fact be 
described as an integral of the halo bias weighted by the halo mass function and by a prescription on how to populate 
halos with galaxies, the Halo Occupation Distribution (HOD). We will return to this issue in Sec. IV C where we will 
discuss the detectability of non-Gaussian corrections to the galaxy power spectrum at large-scales. 



A. Mass function and scale- independent corrections to bias 

The scale-independent contribution can be derived from the halo mass function n(m) with the usual relation [71] 

d In n(m 



b si (m) = 1 + 



38, 



(18) 

<5i=0 



under the assumption of the universality of the mass function and Markovianity in the excursion set derivation, with 
8i representing the large-scale component of the matter fluctuations. In the case of non-Gaussian initial conditions, 
the mass function receives a correction that leads in turn to a scale-independent non-Gaussian correction. Assuming 
the non-Gaussian mass function to be described as nNG(m) = nc{m) RNa{ m , /nl) [73] we can write 

b si (m) = b G (to) + Ab sl (m) , (19) 

where the Gaussian component be is obtained from Eq. (18) using the Gaussian mass function na(m), while the 
non-Gaussian correction is given by [18] 



a 7 / f x d\nR NG {m, / NL ) 
A6 si (m,/ NL ) = — 



88, 



(20) 

a,=o 



We will assume the Sheth & Tormen [74] expression for the Gaussian mass function while for the non-Gaussian 
relative correction Rmg we assume the simple description of [75] based on the Edgeworth expansion of the non- 
Gaussian matter probability distribution function in the Press-Schechter framework [76]. At linear order in /nl this 
is given by 



R NG (m, /nl) = 1 + 7 , x (x 2 - 3) s 3 (x) - - [x - - + G(/ NL ) , (21) 



1 ( 1\ ds 3 (x) 2 

x) rflnT + ° (/NL ' 

where the variable x is defined as x(m) = S c /a m with 8 C = 1.686 being the linear threshold for spherical collapse and 
a m the r.m.s. of the matter perturbations smoothed on the radius R = [3m/(47rp)] 1 / 3 , with p denoting the mean 
matter density. In addition, s 3 (x) represents the reduced skewness of the initial matter density field, Eq. (14). To 
improve the agreement between the measurements of the mass function correction in numerical simulations with the 
prediction of Eq. (21), a scaling parameter q defined by Rng{x) — > Rng(q x ) has been considered [30, 77, 78]. We 
assume here q = 0.91 as derived in [41] from the simulations of [18] (but see also [79], for a similar correction in the 
context of Gaussian initial conditions). 

In the left panel of Fig. 4 we plot the non-Gaussian correction to the halo mass function at z = 0. Both the skewness 
and consequently (in this description) the mass function, at fixed /nl, do not depend strongly on is, at least for v < 1. 
A mild dependence can be noticed as v approaches the limiting value of 1.5. On the right panel of Fig. 4 we show the 
relative scale- independent bias corrections given by Absi/bc as a function of b G for the same models and /nl = 100. 
Such corrections are typically negative and below one percent for moderate values of bo- 
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FIG. 5: Left panel: non-Gaussian correction to the halo mass function at z = 0. Different models are denoted as in the left 
panel. Right panel: Relative, scale-independent bias correction, Abi^i/bi^G, as a function of linear bias 6i,g- 



B. Scale-dependent corrections to bias 

For the scale-dependent correction A&i sc ; to the linear halo bias due to non-Gaussian initial conditions we assume 
the following expression from [29] 



AMm;fc,/ NL )= 2M(fc) ^T- + Mpj°» 



(22) 



where 



I 21 (k,m) = 1 ^ ) J d 3 qM(q)M(\q-k\)W R (q)W R (\q-k\)B„(q,\q-k\ 1 k), (23) 

and with W R (k) a top-hat filter function of radius R corresponding to the halo mass m as m = (4ir /3)pR 3 . This 
description of the linear bias correction is valid under the assumptions of Markovianity and universality of the mass 
function which we adopt here for simplicity. 

For local non-Gaussianity it is easy to show that in the large-scale limit 



I 21 (k,m) k ~° 4/ NL a^, (24) 
so that the second term in Eq. (22) vanishes and the first term gives 

Ab sd (m; k, / NL ) ~ 2 /nl S c ~ 1 , (25) 

that is the usual expression first proposed in [9, 10] with M(k) ~ k 2 for small k. The first term of Eq. (22) however 
generalizes to any non-Gaussian model as in [23], while the second term accounts for additional corrections studied in 
[27, 28]. In addition, the full result of [29] accounts as well for non-Markovian effects and departures from universality 
of the mass function. Such effects will have to be taken into account in a proper comparison with numerical simulations 
results but we can neglect them here as they will not affect our results. 

In the case of the QSF model as described by the template for the curvature bispectrum in Eq. (4), from the 
squeezed limit, Eq. (5), we can derive the large-scale approximation of the integral l2i(k, m), which is given, for v ^ 0, 
by 

t (u \ k ^° f 18\/3~r(i/) ,3/2-!/ v roci-\ 
I 21 (k,m) ~ — /nl 2y -3/ 2 ^(8/27) k S 3/2-„ , (26) 

where, following the notation of [29], 

S„(i?)= / d?qP{q)W 2 R {q)q- n . (27) 
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In the case v = 0, we have 



t n \ £ 36\/3 , 3/, 

/2l(fc ' m) " /NL 23/^y (8/27) fc 



7 + 



ln£) S 3/2 _,-|A^|P(9)^(?) 



(28) 



Since at large scales M(fc) ~ fc 2 , the expected scale-dependent correction to halo bias presents in this case the 
behavior 

AM*) ~ k- 1/2 -" • (29) 

Over the allowed range of values of u, < v < 3/2, such corrections will therefore interpolate between Ab s d(k) <~ k^ 1 / 2 
and A6 s( ;(fc) ~ fc~ 2 , the latter corresponding to the effect of local NG. 

In Fig. 6 we show the quantity l2i(m, k), Eq. (23), as a function of k evaluated for different halo masses and divided 
by its asymptotic value, Eq.s (26) and (28) (see also Fig. 2 and 3 in [29]). We notice that different masses correspond 
to specific dependence on scale for intermediate values of k. In the case of the QSF model for v — 1.5 this aspect is 
more evident than in the otherwise similar local model. It is evident that the asymptotic approximation is able to 
describe the effect only at the very largest scales, while the full evaluation of the integral /2i( m , k) is required already 
for wavenumbers above k ~ 0.01 /iMpc -1 . 

In Fig. 7 we finally show, with thick curves, the relative scale-dependent correction Ab s d{k) /be with b s d given by 
Eq. (22) as a function of k for two representative halo masses m = 10 125 and 10 13 5 h^ 1 M©. Thin curves correspond 
instead to sole contribution due to the first term on the r.h.s. of Eq. (22) while light-colored curves denotes negative 
values. It is evident, in the first place that the case v — 1.5 is indeed very close to the local NG case, up to the 
mentioned numerical factor, both in term of scale-dependence and as a function of mass. As we consider lower values 
of v the dependence on mass, and, in particular, the contribution on the second term on the r.h.s. of Eq. (22) becomes 
relevant, particularly for the low mass examples. In fact, already for v = 1, while in the high mass case such additional 
corrections amounts to a overall decrease in the amplitude of the scale-dependent correction, for the low mass case, 
they induce a change in the sign of the correction itself, albeit keeping a similar absolute amplitude. We should 
remark that this sort of behavior has never been tested in simulations for this specific model, and a dedicated study is 
clearly required. However, as already mentioned, comparisons with different non-local models (e.g. [27, 29]) indicate 
that Eq. (22) does describe correctly the mass dependence of the halo bias corrections to the extent allowed by the 
error on current numerical results. If this description will find further confirmation, possibly for the specific case of 
the QSF model, the results shown of Fig. 7 would be particularly interesting in terms of observational constraints. 
One can imagine, in fact, the possibility of combining the measurements of the power spectrum for different galaxy 
populations where we expect a correction different in sign, allowing in principle for a significant reduction of the 
degeneracy between the non-Gaussian parameters and the linear Gaussian bias (see e.g. [43, 49] for recent works in 
this direction). The other panels of Fig. 7 referring to the QSF model, show that a similar situation is recovered for 
even smaller values of v, as for instance v = 0.5 corresponding to a scale-dependent correction Ab s d ~ 1/k. 



C. Fisher matrix 



We perform a Fisher matrix analysis in order to determine the expected, simultaneous constraints on both the 
/nl and v parameters from future galaxy surveys. Similar forecasts, in terms of the single /nl parameter, have 
been performed, for the local model alone in [9, 42, 46-48] for upcoming galaxy and clusters surveys. The cases 
of equilateral and orthogonal non-Gaussianity, in addition to the local model, have been studied by [50] and [51], 
focusing respectively on the EUCLID survey [80] (combining weak lensing with photometric and spectroscopic data) 
and the eRosita, X-ray cluster survey [81]. In these works, the effects of non-local models include scale-dependent 
corrections to halo bias, described by expressions analogous to Eq. (22). Running non-Gaussianities described by an 
additional running parameter ting, such as those for both the local and equilateral shapes [53, 54], also depend on 
two parameters. Forecasts for these models have been considered in [50, 52, 75]. Some examples of these running non- 
Gaussianities may also alter the momentum dependence in the squeezed limit of the bispectrum. But the difference 
between this case and the QSF inflation is clear. The former is caused by the running of non-Gaussianity and the latter 
by the shape of non-Gaussianity. For QSF inflation, even if 3/2 — v is not much less than one and so the deviation 
from the local shape is significant, the bispectrum can be still scale-invariant. But for running Gaussianities, such a 
momentum dependence corresponds to a case with very strong overall scale-dependence. 

Our goal is to provide an estimate of the possibility to constrain the parameter v assuming a positive detection 
of non-Gaussianity by Planck. For this purpose we consider two large-volume surveys comparable in size and red- 
shfit range to EUCLID and LSST. Such surveys are already expected to provide constraints to f^£' [46, 50] from 
measurements of the galaxy power spectrum comparable to those expected by Planck [82-84]. Being this a first, 
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FIG. 6: Comparison of the full calculation of the function I21, Eq. (23), with its asymptotic value as k — > for the QSF model 
with v = 1.5 (top left panel), v = 1 (top right panel), v = 0.5 (central left panel) and u — (central right panel). The lower panels 
show the same quantities for the local (left) and equilateral (right) models. Dotted, dashed and continuous curves correspond 
respectively to halo masses m = 10 12 , 10 13 and 10 14 hT 1 Mq. Notice that the plotted ratio is independent of redshfit. The 



noisy results for the v ■■ 
terms on the r.h.s. 



case are due to the vanishing of the asymptotic expression, Eq. (28) due to the canceling of the two 



indicative assessment, we will assume a simplified description, characterized simply by the field of view, the redshift 
range (i.e. the volume) and the expected linear (Gaussian) bias b(z) as a function of redshift. Since we expect most of 
the signal to come from the largest scale probed, we do not expect a significant impact of shot-noise (see, for instance, 
[46, 50]). On the other hand, shot-noise should not affect, by design, the BAO analysis, the primary target of these 
missions. We will nevertheless include a shot-noise contribution to the power spectrum variance as detailed below. 
The photometric or spectroscopic nature of redshfit observations is assumed to play a negligible role as, again, the 
determination of large-scale power does not require high precision in the radial galaxy positions. 

For our analysis we assume a ( "EUCLID-like" ) survey, denoted as VI, with a 20,000 deg 2 field of view and a redshift 
range of 0.4 < z < 2, and a galaxy population with fiducial (Gaussian) linear bias parameter given by bc(z) = \fl + z 
[85]. This would allow a comparison with the results of [50] for the local and equilateral models. As opposed to [50], 
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FIG. 7: Relative scale-dependent correction Ab s d(k) jba (thick dashed and thick continuous curves) as a function of k for two 
representative halo masses m = 10 12 ' 5 Mq (dashed ones) and m = 10 13 ' 5 Mq (continuous ones), among which the light- 
colored ones denote negative values. The choice of the models shown is the same as Fig. 6. Thin dashed and thin continuous 
curves correspond to sole contribution due to the first term on the r.h.s. of Eq. (22). 



which considers 12 equally populated redshift distributions, we will simply compute the Fisher matrix information in 
redshift bins of size Az = 0.1. As an exercise we extend this galaxy population, described by the same bias evolution 
bc{z), to a larger ("LSST-like") volume, denoted as V2, given by a field of view of 30,000 deg 2 and redshift range 
0.3 < z < 3.8. This corresponds to the second example of [46], although with a lower (therefore more conservative) 
fiducial value for the linear bias. 

We consider a two-dimensional Fisher matrix F a ,b for the parameters p a = {/nl, v} obtained as a sum over all N z 
redshift bins of the three-dimensional matrix for the parameters p a = {/nl ; v t ^g}: marginalized over the linear 
bias parameter be- Since the scale-dependent correction Ab s d depends linearly on the product /nl Q>g ~ 1)> we can 
expect a significant degeneracy with /nl and a marginalization over the value of b(j should be taken into account. 
Such marginalization is performed in [50] but not in [46], where, on the other hand, only the very largest scales are 
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considered for the analysis of the local model alone. The matrices are defined respectively as 



(30) 



and 



*(z 3 ) 



J~ rv 



E 



dPgik^ZijdPgikj 



1 



.(2 3 ) 



Pa 



d Pb AP*(kj,Zi) 



(31) 



where the sum runs over the available wavenumbers from the fundamental frequency of the rcdshfit bin k n 
2tt/V 1 / 3 (z), V(z) being the bin volume to /c max , in steps of k m i n . 
The expression for the galaxy power spectrum is given by 

P g {k,z) = b 2 (z)P(k,z) 

= [b 2 G (z) + 2 b G {z) (Ab si (z) + Ab sd (k, z))] P(k, z) + 0(f* L ) 



(32) 



where we include only linear corrections in / NL . It should be noted that primordial non-Gaussianity modifies both the 
amplitude and the form of the galaxy power spectrum and a large range in k helps reducing degeneracies between non- 
Gaussian parameters and bias. Following [50], we choose & max = ir/(2 R) with the scale R obtained from the redshift- 
dependent equality u 2 R {z) = 0.48 constant. At z = this choice provides, for our cosmology, fc max = 0.15 /iMpc -1 , 
as in [50] and it implies the range of values 0.2 < fc max < 0.5 for 0.5 < z < 2. We will later discuss how our results 
depend on this choice, comparing them with those obtained from a more conservative assumption corresponding 
to fc max = 0.075 /iMpc -1 at z = 0. For the matter power spectrum P(k,z) we assume the linear expression and 
consequently ignore nonlinear corrections at small scales. While on one side we do expect such corrections to be 
there for k close to fc max , the final results do not strongly depend on this assumption since we neglect them as well 
in the power spectrum variance AP 2 . In fact, in our approximation, their inclusion would only reduce the impact 
of shot-noise on AP 2 The range in wavenumbers is, on the other hand, well within reach of accurate predictions in 
nonlinear perturbation theory (see e.g. [86]). Notice that, for simplicity, we also ignore possible corrections due to 
nonlinear bias. This is in part justified by the small value of the quadratic bias parameter corresponding to the linear 
one b G (z) (typically \b 2 /b\ < 0.35) that can be derived by the expressions of [87]. 
The power spectrum variance is given by (see e.g. [52]) 



AP 2 (k,z) 



2irk 2 



1 + 



4Ab NG (k,z)P(k,z) 



Ptot{k,z) 



+ 0(/nl) 



(33) 



which accounts for non-Gaussian corrections to the bias and where P to t is the Gaussian galaxy power spectrum 
including shot-noise, 



P to t(k,z) =P g (k,z) + 



(27r) 3 n(z,) ' 



(34) 



Since we are describing the galaxy population simply in terms of the value of the linear bias we derive the mean 
number density n(zi) in the redshift bin Zi making the quite drastic but simplifying assumption that each dark matter 
halo of mass above a certain minimal mass M min {zi) contains a single galaxy. The value of M min (zi) is determined 
imposing the relation 



b G ( Zl ) = VT+Y, = ( 

J A 



n G (m, z^ b G ,h(m, z t )mdm, 



(35) 



where we assume the Sheth-Tormen [74] expressions for the Gaussian halo mass function n G (m,z) and linear halo 
bias b Gy h{in, z). The galaxy density is then obtained as 



n(zi) = / n(m,Zi)dm. 

JM min {zi) 



(36) 



We notice that since the halo bias b Gi h(m,Zi) weights more massive halos with respect to smaller mass halos, this 
procedure tends to overestimate the value of M m i n and therefore underestimates the value of the density when 
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TABLE I: Marginalized 1-a errors on the parameters /nl and v expected for our two example surveys VI and V2 for a choice 
of fiducial values. For a comparison with other works, we include the constraints expected for the local and equilateral models. 



compared to the same calculation performed assuming the proper Halo Occupation Distribution (HOD) for the 
galaxy sample. In addition, the definition of the minimal mass M min allows us, perhaps improperly, to define a mean 
or "characteristic mass" for the halo population of the bin at Zi given by 

fhi = / n(m,Zi)mdm. (37) 

JM min {zi) 

We assume the mean mass fh(zi) to evaluate the integral l2i(k,m,z), Eq. (23), and its derivatives for each bin. 
Notice that in the marginalization over the bias, i.e. in the derivative dP g /dba, we not only consider the explicit 
dependence on the parameter be but derive as well all mass-dependent quantities like a m or the integrals as li\ 
as d/db G — [db G ^ h (m, z)/dm]~ 1 {d/dm). Obviously, the correct procedure would have involved integrating all bias 
corrections over the proper range of halo masses and the proper HOD. However, such a drastic solution is still more 
conservative than neglecting altogether the marginalization over the mass-dependence of such corrections as done 
in previous similar analysis in the literature. Ref. [50] correctly points out that while the mass-dependence of such 
corrections is relevant in the case of the equilateral model, predictions could not be properly tested in simulations yet 
for this specific model (see e.g. [29]) and they restrict their expressions for Ab s d to its asymptotic value at small k. 
However, we assume the full expression in Eq. (22) to be valid, noticing that such issues are less important for our 
model when v > 0.5 and that after marginalization over bias (and mass), our results are, as we will see, consistent 
with those of [50] for both the local and equilateral models. 

Finally, we do not consider a full marginalization over cosmological parameters. Forecasted constraints on the local 
/nl' parameter, from measurements of the galaxy power spectrum marginalized over the cosmology (with priors from 
the Planck CMB power spectrum) have been studied in [46] finding an increase in the error A/#£- of about 30% for 
a EUCLID-like experiment. Ref. [50] compares instead errors on different NG models marginalized on cosmology 
with and without Planck priors. They find, as one can expect, a particularly strong degeneracy of the equilateral 
/nl parameter with cosmological ones, due to the lack of a scale-dependence in the halo bias correction. We might 
therefore expect an larger impact of the uncertainty on cosmological parameters when low value of v are considered. 

The marginalized 1-a errors on the /nl and v parameters obtained for the VI and V2 geometries assuming both 
J: mal = 0.15 /iMpc" 1 and fc max = 0.075 h Mpc -1 at z = are shown in Table I. In the first place we notice that 
uncertainties obtained for the local and equilateral parameters are consistent with the results of [50]. We find in fact 
lower errors respectively by 30% and 50% with respect to their uncertainties when marginalized over cosmological 
parameters with Planck priors. In the case of the QSF model we observe the expected improvement of the constraints 
when larger values of both /nl and v are assumed as fiducial values. The dependence on the fiducial values of /nl 
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FIG. 8: 1- and 2-a uncertainty contours corresponding to the determination the parameters /nl and v for three different choices 
of the fiducial values: /nl = 50 (left column), 100 (central column) and 150 (right column) and v — 1.5 (top row), 1.0 (central 
row) and 0.5 (bottom row). Continuous blue curves corresponds to the VI geometry, dashed, red ones to the V2 one. Both 
cases assume a limiting fc ma x(z) such that fc max (0) = 0.15 ftMpc -1 . Caution should be taken when the elliptical contours are 
large (see discussions in Sec. VB.) 



and v of the 1-cr uncertainty Ah* can be very roughly described over the range 0.5 < v < 1.5, for the VI geometry, as 

0.1 100 

A^~^— , (38) 
v 6 /nl 

with the V2 case corresponding to an uncertainty about a factor of two smaller. The important results is that we can 
indeed expect to able to distinguish values of v ~ 1 from the local limit v = 1.5 for sufficiently large values of /nlj 
say /nl > 100, already with a survey corresponding to our VI example. Even lower values of /nl would be sufficient 
for a larger survey as V2. The improvements in the constraints in the V2 example with respect to VI are mainly due 
to the larger volume, although the larger values of the fiducial Gaussian bias and the higher redshift do play a non 
negligible role. 

These considerations are particularly evident in Fig. 8, where we show the 1- and 2-a uncertainty contours for /nl 
and v corresponding to the VI (blue, continuous curves) and V2 (red, dashed curves) examples, assuming the set of 
fiducial values given by /nl = 50, 100 and 150 and v = 1.5, 1.0 and 0.5. Notice that we choose to keep the same 
range over the variable v for the case of fiducial v = 1.5 as for the other values considered. While such choice makes 
it hard to distinguish the different curves for large fiducial value of /nl and v, it allows, on the other hand, an easier 
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comparison of the results as the fiducial v is varied. The complete results of the Fisher analysis are given, in any 
event, in Table I. 

The degeneracy between /nl and v can be easily understood for values of v close to the local limit v = 1.5. In 
this case, in fact, since the dominant effect is given by the scale-dependent correction Ab s d and a lower value of 
v can be compensated by an higher /nl. When v is close to one, the degeneracy is reduced since the large-scale 
scale-dependent corrections are balanced by almost scale-independent corrections at small scales. Notice that the 
bias corrections shown in Fig. 7 are evaluated at z = 1, an at this redshfit the relevant range of scale is given by 
0.003 hMpc^ 1 < k < 0.28 ftMpc" 1 . At lower values of v the error Ai/ increases significantly, as the scale-dependent 
corrections are now very mild and there is a large degeneracy between all the parameters in each redshift bin. The 
combination of different bins still allows for a relatively low uncertainty on /nl- 

Fig. 9 show the same results as Fig. 8 but for a limiting fc max (z) such that fc max (0) = 0.075 h Mpc -1 . We notice, 
in the first place, a significant worsening of the determination of /nl already at v = 1. This is due to the lower 
control over the form of the power spectrum: the constrains are now mostly based on its amplitude at large-scales. 
Remarkably the overall increase in the uncertainty on v is only of about a factor of two for all considered fiducial 
values. The different choice of the range of scales assumed for the analysis results as well in a different degeneracy 
between /nl and v. This is evident for instance for the fiducial values /nl = 100 and v = 1. In the conservative 
case of A: max = 0.075 /iMpc -1 at z — 0, in fact lower values of v are degenerate with higher values of /nl as expected 
from the scale-dependence alone. Including smaller scales the degeneracy is greatly reduced and turns slightly in the 
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opposite direction. We should stress the fact that the possibility to properly include in the analysis scales in the mildly 
nonlinear regime relies on accurate predictions for both the matter power spectrum and bias. Recent developments in 
the description of matter nonlinearities in Perturbation Theory (see e.g. [88-90]) can justify even our more optimistic 
choice of fc max = 0.15 h Mpc -1 at z = 0, also in the case of non-Gaussian initial conditions [91, 92]. A similarly 
accurate description of halo and galaxy bias in this regime is a more challenging task and will necessarily require 
further work (see e.g. [93, 94] for recent results). 

We finally remark that a similar analysis can be extended to the galaxy bispectrum. As shown in [21, 22, 38-40] 
a scale-dependent correction as the one considered here for linear halo bias is present as well for terms induced by 
nonlinear, quadratic bias and relevant for the galaxy bispectrum at large scales. A relatively simple model for the halo 
bispectrum in the presence of local non-Gaussian initial conditions has been recently tested in numerical simulations 
by [41] where the addition of scale-dependent bias corrections are shown to improve significantly over previous studies 
of the galaxy bispectrum [12, 13]. 



V. CMB CORRELATIONS AND JOINT CONSTRAINTS 



We now turn to consider the constraints on Quasi-Single Field inflation from the bispectrum of the cosmic microwave 
sky, notably with the aim of making forecasts relevant for the Planck data. In principle, given the linearity of the 
CMB transfer functions, this is more straightforward than for the matter perturbations. In practice, however, the 
task is very computationally intensive so we shall make some simplifying approximations which will yield constraints 
and forecasts of reasonable precision. 

Our aim is to compare a theoretical prediction for the CMB bispectrum with that obtained from observations, such 
as WMAP or Planck full sky maps. The CMB bispectrum is the three-point correlator of the harmonic coefficients 
a lm describing the map, B m l ^ 2Tri3 = a; 1 m 1 0'i2m2 a i3m 3 - Assuming the bispectrum has been created by a statistically 
isotropic process, we can restrict our attention to the angle-averaged bispectrum Bi 1{ 



2'3 > 



Bl t l 2 l 3 ~~^^hl 1 l 2 l 3 Gmim 2 rn3 a hm 1 CI'l2rn2 a hm 3 , (39) 

where hi 1 i 2 i 3 is a geometrical factor which enforces a multipole triangle condition, 



_ ./ (2 il + 1 )(2 ii + l)(2 fa + l) h h h , (40) 



4tt I 



and Q m \ml m is the Gaunt integral, 



nhhh =h I h l 2 h \ 
1 mi m 2 m 3 



with the usual Wigner-3j symbol. It is generally more straightforward to work with the reduced bispectrum 6; 1 ; 2 z 3 = 
i B^i^ with the geometrical factors removed. 
To find the reduced CMB bispectrum 6i 1 ; 2 ; 3 induced by a given primordial bispectrum Bo(ki,k 2 ,k 3 ) 1 we use the 
CMB transfer functions Aj(fc) to project forward as 

b hhh = (J^J Jdxdk 1 dk 2 dk 3 (xkik 2 k 3 ) 2 B {ki, k 2 , h) 

A h {h)A h (k 2 )A h (fc 3 ) 3h {kix)ji 2 {k 2 x)j h (k 3 x) . (42) 

If the original bispectrum B^{ki,k 2l k 3 ) is separable, the complicated 4D integral (42) also separates and becomes 
much more tractable. While the local model and the usual equilateral ansatz are separable, the quasi-local shape, 
Eq. (4), which interpolates between them is is not separable. Nevertheless, for a given QSF parameter we can still 
calculate 6; 1 / 2 i 3 by using a separable eigenmode expansion, but we will not describe the modal methodology here in 
detail (see [95]). 
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1.0 




v 



FIG. 10: Left panel: Estimate of the CMB correlation of the quasi-single field shape as a function of the parameter v with 
the local and equilateral shapes using the weighted shape correlator Eq. (46). Note that the v = 1.5 model is nearly perfectly 
correlated with the local shape, whereas an 88% correlation is achieved in the equilateral limit v — > 0. Right panel: Estimate 
of the CMB correlations of the QSF models with themselves C(B,B') as a function of the parameters v, v' . Note the large 
plateau of strong correlations for models with < v < 0.75 and < v' < 0.75, which makes it difficult distinguish QSF models 
with small v. 



A. CMB bispectrum correlator 



In order to determine if a given theoretical bispectrum is present in the observational data, we employ the approx- 
imate estimator 



S = 



1 



frfo c h c l2 c l3 



(43) 



with appropriate experimental effects incorporated - beam, noise and mask. The estimator essentially performs a 
least squares fit between theory and data with the ratio above yielding the signal-to-noise. We can, in principle, use 
the separable modal methodology in ref. [95] to find constraints on quasi-single field inflation. The reconstructed 
bispectrum coefficients extracted from the WMAP7 data [95] can be used to compare directly with the expansion 
coefficients predicted theoretically from Eq. (42). 

Whether or not two values of the quasi-single field bispectrum parameter v can be distinguished by a given data 
set can be determined by a Fisher matrix analysis, essentially the cross-correlator between the two CMB bispectra 



B, 



B' 



(44) 



where the normalization is defined by 



N 



E 



\ ^ Ci x Ci 2 Ci z \^ *j* t Ci x Ci 2 Ci 



E 



B> 



(45) 



However, this approach is extremely computationally demanding as we must calculate the full bispectrum for each 
value of v before we can make any comparison. This has been achieved already for local and equilateral asymptotes 
but we would like a simpler method for estimating the interpolants. As shown in [96], a fairly accurate measure of 
the Fisher matrix, Eq. (44), can be obtained from the shape correlator, 



C(S,S') = 



F(S,S') 
y/F(S,S)F(S',S') 



(46) 



where 



F(S,S')= / S(k 1 ,k 2 ,k 3 )S'(k 1 ,k 2 ,k :i )u J (k 1 ,k2,k 3 )dV k , 



(47) 
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with the bispectrum 'shape' S(k\, k 2 ,k 3 ) defined from the rescaled primordial bispectrum 

S(h, k 2 , k 3 ) = (fcifc 2 fc 3 ) 2 B*(fci, k 2 , h) . (48) 

Here, in order to replicate the scaling in the CMB bispectrum correlator, Eq. (44), we adopt the weight function in 
Eq. (47) as 

It is demonstrated in ref. [96] that Eq. (46) yields a good phenomenological approximation to Eq. (44) for a wide variety 
of shapes (all classes of scale-invariant models). For example, the shape cross-correlator C(S,S') for the local and 
equilateral models, produced only a 5% underestimate of the true CMB correlator C(B, B') at Planck resolution. This 
essentially sets an upper limit on the accuracy of the weighted shape correlator for estimating the CMB correlations 
between the local, equilateral and QSF models. It is sufficient precision for making Planck forecasts in the present 
context. 

Estimates of CMB correlations for the QSF models using the shape correlator of Eq. (46) are illustrated in Fig. 10 
(left panel). It is clear from Fig. 10 that the QSF model with v = 1.5 is very well approximated by the local shape, 
while smaller values across a broad range v < 0.75 are close to equilateral, though they never attain more than a 
90% correlation. The two-dimensional cross-correlation for the QSF models with themselves (i.e. as functions of v, 
v') from Fig. 10 (right panel) indicate that it will be very difficult to distinguish between models with v, v' < 0.75. 
However, there are significant differences between models with small v < 0.75 and large v' > 1.0. 



B. CMB Fisher matrix and likelihood analysis 

In order to forecast the implications of the Planck experiment for QSF models, we have performed a two-dimensional 
CMB Fisher matrix analysis analogous to that for the large-scale structure in Section IV C, as well as a likelihood 
analysis. Here we have defined a CMB Fisher matrix T a ,b f° r the parameters p a = {./nl, ^} as in Eq. (31) using 
knowledge of the approximate correlator, Eq. (46), for the QSF shape functions. We do not marginalize over cos- 
mological parameters for which we believe there is only a weak dependence for the CMB bispectrum in any case. 
Assuming that Planck achieves a target 1-a variance A/^l' = 5, in Figure 11 we have plotted l-a and 2-er uncertainty 
contours for /nl and v for the same fiducial values as previously: /nl = 50, 100 and 150 and v — 1.5, 1.0 and 0.5. 
These CMB Fisher matrix ellipses (continuous curves) have degeneracy directions which are complementary to the 
large-scale structure analysis, and so there is value in a joint analysis as we shall discuss. These results show that 
v i=a 1.0 models will be distinguishable from both equilateral or local models at better than the 2-a level for /nl > 100 
(which is approximately equivalent when normalised relative to the local model to /^l' > 30). Nevertheless the results 
also show the deficiency of the Fisher matrix ellipses for models near v 1 v' < 0.75 where the correlation functions are 
very flat or degenerate. The large elliptical contours obtained from the local derivatives around v = 0.5 appear to 
indicate that these models cannot be distinguished from the local models, which is apparently not consistent with the 
v > 1.0 results at large /nl- 

The shortcomings of this simple Fisher matrix analysis motivated an improved likelihood analysis based on the 
approximate correlator, Eq. (46). This was achieved in several steps. First, the Fisher matrix T for n=300 values 
of v (denote these discrete values Vj) was calculated on a uniform n 2 grid using Eq. (47). All matrix elements 
were divided by the local-local correlation result, so that the diagonal elements of this matrix yield the QSF model 
normalization at a given Vj relative to the local bispectrum. Of course, when normalised with Eq. (46) then T yields 
the correlation matrix C illustrated in Fig. 10 (right panel). Secondly, we diagonalize the Fisher matrix to extract the 
optimal uncorrelated basis, that is, representing the matrix as 

T = VVV T (50) 

where V contains the unit eigenvectors Vi and I? is a diagonal matrix constructed from an ordered list of the eigenvalues 
Xi. PCA analysis of T for these QSF models reveals that just two eigenvalues are dominant, with eigenvectors 
representing the local shape and an equilateral-like shape. A third eigenvector makes a relatively small additional 
contribution for a more accurate analysis. 

We next perform Monte Carlo simulations to populate data distributions about a specific {/nl, ^j} model which 
are consistent with the Planck variance A/^l' = 5. To achieve this, we need to know how much of each eigenvector 
Vi is required to build a particular correlation curve, corresponding to the given i/j value. The corresponding data 
'shape' vector for the Vj model can be found as the j th row of the matrix 



(51) 
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FIG. 11: Expected joint constraints on the two parameters /nl and v from measurements of temperature bispectrum from 
Planck-like CMB observations. The shaded areas of decreasing colour density corresponds to the 1-, 2- and 3-cr constraints 
from the forecasted likelihood function. For comparison, the Fisher matrix results corresponding to 1- and 2-a constraints are 
show by superimposed ellipses. These are qualitatively in agreement and identify the correct degeneracy directions, but become 
inaccurate in degenerate regions for small v. 



which is automatically unit normalised. We then multiply by the required /nl (again normalised to the local model) 
and the appropriate noise vector Nj is added, which consists of n Gaussian random variables of variance 5, that is, 
the data vector becomes 

* = /nl ^ + Ni . (52) 

We can do this essentially because we are working in an orthonormal 'observation' space in which the probability dis- 
tribution is a multivariate Gaussian (this understanding of the independent variance in orthogonal modes is exploited, 
for example, in modal CMB map simulations described in ref. [95]). We must then map back to the {/nl, v) space 
in which /nl (normalised to the local model) for each given original value of Vj is given by /nl = diCe^ 3 and the 
best fit v for this data is the one that minimises the error X^(^i — /nlo^) 2 - This process is repeated many times 



22 




FIG. 12: Left panel: Minimal values of /nl required to separate the QSF shape from either the spurious local or spurious 
equilateral signal it would generate at different significance levels. The black continuous curve shows the value of /nl for which 
either the spurious local or spurious equilateral signal would violate the WMAP7 2-<r bounds on these parameters. Right panel: 
Expected joint constraints on the two parameters /nl and v from measurements of temperature bispectrum from Planck-like 
CMB observations, assuming the fiducial value /nl = 0. The shaded area of decreasing colour density correspond to the 
1-, 2- and 3-<j constraints from the forecasted likelihood function. Note that volume factors in the transformation from the 
orthonormal 'observation' space back to the {/nl, ^}-space yields a non-uniform distribution for v which disfavors values of 

to build up a smooth likelihood function about the given model. This mapping from the orthonormal 'observation' 
space back to {/nl, v\ space does not necessarily recover a peak in the distribution around the original model values 
because the Jacobian for the transformation rescales the volume factor. The non-trivial weighting in the {/nl, v\ 
'theory' space introduces artifacts like those shown in Fig. 12 (right panel) where v does not even have a uniform 
probability distribution about the null /nl = model. This indicates that it can be useful to choose parametrisations 
for theoretical models which match, as closely as possible, the signal-to-noise weighting mandated by the observational 
data. 

The results of the CMB likelihood analysis are illustrated in Fig. 11 with l-er, 2-a and 3-er contour shading. A direct 
comparison to the Fisher matrix ellipses is made for the same input model values as above. While the results are 
in qualitative agreement and identify the same degeneracy directions, they also illustrate the limitations of a Fisher 
matrix analysis which requires a 'global' extrapolation into regions of the parameter space where the distributions are 
either rapidly changing or degenerate, neither of which can be adequately captured by local derivatives. Additional 
features appear in the likelihood analysis, particularly note the two peaks in the nearly degenerate v = 0.5 case. For 
small v, the likelihood contours are cut off before reaching the local model (v = 1.5), so there are parameter values 
for the QSF model which can be distinguished from both local and equilateral shapes for f^ L F > 100. 

In Fig. 12 (left panel) we illustrate the expected /nl signal required to distinguish QSF from the local and equilateral 
models at various significance levels. We also plot the present WMAP7 constraints on any QSF model from the union 
of local model and equilateral model estimation results [7] (assuming the correlations shown in Fig. 10). As can be 
seen, it will be difficult to distinguish equilateral from QSF models with v < 1.0 where the correlation function C is 
very flat; a large measurement /nl > 200 is required and this is excluded already by WMAP7. Similarly, it will not be 
possible to distinguish the local model from QSF models in the narrow region above v > 1.3 because the necessary /nl 
is also excluded. This leaves a narrow window of opportunity for Planck data to identify QSF models for 1.0 < v < 1.3 
assuming a large /nl > 120 measurement (approximately equivalent to a local /j^' > 40). Using Planck polarization 
data with a reduced variance A/nl ~ 3, the forecast window widens considerably to 0.5 < v < 1.4 and offers a much 
better prospect of differentiating between the QSF and other models at greater than a 2a significance. 

C. Combined CMB and LSS forecasts 

In this section we consider the combination of our results for the LSS Fisher matrix analysis with those from the 
approximate CMB analysis. This combination is particularly interesting because, as mentioned above, the CMB 
bispectrum and the galaxy power spectrum observations weight differently distinct triangular configurations of the 
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FIG. 13: 1- and 2-a CMB contours (black, dot-dashed curves), LSS contours corresponding to the VI survey, assuming fc max = 
0.15 h Mpc -1 at z — (blue, continuous, thin curves) and the contours corresponding to the joint CMB and LSS analysis (blue, 
continuous, thick curves). As in the previous plots we consider the set of fiducial values given by all combinations of /nl = 50, 
100 and 150 and v = 0.5, 1 and 1.5. 



initial bispectrum. For this reason we expect that a joint forecast should reduce the specific degeneracies between the 
two parameters /nl and v already highlighted in the previous sections. 

Since we did not consider a full likelihood analysis for the LSS forecasts, we will limit ourselves to comparing the 
Fisher matrix results. Fig. 13 shows the (1- and 2-a) CMB contours, the LSS contours corresponding to the VI survey, 
assuming & max = 0.15ft.Mpc _1 at z = 0, and the contours corresponding to the joint CMB and LSS analysis. The 
complementarity of the two observables is evident, for instance, in the case given by the fiducial values /nl = 100 
and v = 1 (central panel). We remark, again, that the large-scale structure results rely on our ability to accurately 
describe non-Gaussian correction to galaxy bias at mildly nonlinear scales, where the scale-independent component 
of such corrections dominates. This can be seen, for instance, from the particular degeneracy between /nl and v. 
As already mentioned, an accurate theoretical model for galaxy bias at these scales will require further investigations 
and comparison with numerical simulations. 

The same joint results, but assuming the more conservative choice of fc max = 0.075 h Mpc" 1 at z — for the 
galaxy power spectrum analysis, are shown in Fig. 14. In this case the scale-dependent correction to bias represent 
the dominant effect of primordial non-Gaussianity on the LSS observable. Remarkably, for a fiducial v = 1, while the 
uncertainty on /nl from the galaxy power spectrum is very large, the addition of LSS information can significantly 
improve the CMB constraints on v. 



24 




Table II shows the marginalized l-er errors predicted by the CMB Fisher matrix analysis. In addition, also shown 
is the combination of the CMB analysis with the LSS results for the VI and V2 surveys assuming the two choices for 
the range of scales included considered in Section IV. 



VI. CONCLUSIONS 



In this work we have studied the effects of the primordial non-Gaussianity predicted by Quasi-Single Field models of 
inflation on large-scale structure observables with particular attention to the galaxy power spectrum and to the induced 
scale-dependent correction to halo and galaxy bias. In addition we studied the detectability and the correlation with 
other non-Gaussian models in CMB bispectrum measurements. 

Scale-dependent bias corrections have been the subject of several studies in recent years as they provide a remark- 
able test of squeezed configurations of the initial bispectrum. In the particular case of the local model, where the 
primordial curvature bispectrum peaks precisely in the squeezed limit, constraints on the amplitude parameter /nl 
from current measurements of the power spectrum of high-redshift sources are already comparable to CMB results. 
In this context, Quasi-Single Field models represent a veritable case study as they predict a one-parameter family of 
curvature bispectra with variable momentum-dependence in the squeezed limit, resulting in a correction to linear bias 
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TABLE II: Marginalized 1-a errors on the parameters /nl and v expected from the Fisher matrix analysis of the Planck 
temperature bispectrum alone and in combination with the two example LSS surveys VI and V2 for a choice of fiducial values. 

with a scale-dependence interpolating between the one predicted by the local model to practically scale-independence. 

In addition to this interesting phcnomcnological aspect, the determination of the parameter v can also provide direct 
information on the high-energy theory. Supersymmetry naturally determines the range of masses of light scalars during 
inflation. Interestingly, such masses may be directly observable through QSF inflation models, in terms of the powers 
in the peculiar scale-dependence in galaxy bias and/or momentum-dependence in CMB non-Gaussianitics. In this 
paper we investigated how much, in the event of the discovery of large primordial non-Gaussianities, we can refine 
such discovery and tell the difference between QSF predictions from those of the local or equilateral shapes in the 
future LSS and CMB experiments. 

We find that there is an interesting cosmic complementarity between the galaxy halo bias and CMB in the two- 
parameter space, v and /nl ■ The halo bias is relatively more sensitive to v since the scale-dependent correction mostly 
depend on the squeezed-limit behavior of the primordial bispectrum. The CMB is more sensitive to /nl since it has 
more signal weight away from the squeezed limit. By combining both, depending on the value of /nl, & significant 
fraction of the parameter space in v may be distinguished from the local model, although for v very close to 3/2 the 
degeneracy will remain. 

In this work we have concentrated on the primordial bispectra. Nonetheless QSF models also predict large 
trispectra with similar peculiar momentum dependence. Furthermore, at least in the perturbative region, these 
models predict that the size of the trispectra (from the scalar-exchange diagram) is much larger than the size 
of the bispectra squared / NL [2]. In particular, various soft momentum limits of the trispectra have corresponding 
intermediate momentum dependence determined by the isocurvaton mass [97]. These make the observability of 
trispectra more prominent comparing to models with t^p L <~ / NL [98-100]. It remains to be seen, however, whether 
the promise offered by the trispectra can contribute further towards the goal of this paper. 

During the completion of this project we became aware of a similar study in preparation [101]. A rough comparison 
with their results for the galaxy power spectrum analysis indicates a broad agreement. We are grateful to the authors 
for coordinating with us the publication on the ArXiv website. 
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